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This report shows the latest results on the study of the method to determine the angle 4*3 of 
the unitarity triangle using Dalitz plot analysis of D° decay from B — > DK ± process in a model- 
independent way. We concentrate on the case with a limited charm data sample, which will be 
available from the CLEO-c collaboration in the nearest future, with the main goal to find the 
optimal strategy for 03 extraction. 



I. INTRODUCTION 



II. MODEL-INDEPENDENT APPROACH 
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The measurement of the angle 03 (7) of the unitarity 
triangle using Dalitz plot analysis of the D° — > KgTr + ir~ 
decay from B ± — > DK ± process, introduced by Giri et 
al. [H and Belle collaboration 0] and successfully imple- 
mented by BaBar Q and Belle [J], presently offers the 
best constraints on this quantity. However, this tech- 
nique is sensitive to the choice of the model used to de- 
scribe the three-body D° decay. Currently, this uncer- 
tainty is estimated to be ~ 10° and due to large statistical 
error does not affect the precision of 03 measurement. As 
the amount of B factory data increases, though, this un- 
certainty will become a major limitation. Fortunately, a 
model- independent approach exists (see [l|), which uses 
the data of the r-charm factory to obtain the missing 
information about the D° decay amplitude. 

In our previous study of the model-independent Dalitz 
analysis technique [j| we have implemented a procedure 
proposed by Giri et al. involving the division of the Dalitz 
plots into bins, and shown that this procedure allows to 
measure the phase 03 with the statistical precision only 
30-40% worse than in the unbinned model-dependent 
case. We did not attempt to optimize the binning and 
mainly considered a high-statistics limit with an aim to 
estimate the sensitivity of the future super-B factory. 

The data useful for model-independent measurement 
are presently available from the CLEO-c experiment @. 
CLEO-c collected an integrated luminosity of 280 pb _1 
at the ?/>(3770) resonance decaying to DD. By the 
end of CLEO-c operation this statistics will grow up to 
750 fb _1 0. This corresponds to ~ 1000 events where 
D meson in a CP eigenstate decays to Kgir + ir~, and 
twice as much events of t/03770) — > D°D° with both D 
mesons decaying to K°ir + ir~~ . Both of these processes in- 
clude the information useful for a model-independent 03 
measurement. In this paper, we report on studies of the 
model-independent approach with a limited statistics of 
both ^(3770) and B data, using both Dcp — > KgTr + n^ 
and (K°TT + ir~)D(K a 7r + Tr~)E) final states. 



The density of D° — > Kg7r + Tr Dalitz plot is given by 
the absolute value of the amplitude jo squared: 

Pd = PD{m+,m 2 -) = |/d(to+,W-)I 2 (1) 

In the case of no CP- violation in D decay the density of 
the D° decay Pd equals to 



Pd 



\Id\ 2 = PD(m 2 _,m 2 + ). 



(2) 



Then the density of the D decay Dalitz plot from B^ 
DK ± process is expressed as 

Pb± =\fD + r B e^ B± ^f D \ 2 = 

Pd + r 2 B p D + 2 y / p D p D (x ± c + y±s), 



(3) 



where x±,y± include the value of c/>3 and other related 
quantities, the strong phase 8b of the B ± — > DK ± decay, 
and amplitude ratio tb'- 

x±=r B cos(S B ±(f>3); y± = r B sm(S B ± cf> 3 ). (4) 

The functions c and s are the cosine and sine of the strong 
phase difference ASd between the symmetric Dalitz plot 
points: 



c = cos(<5£)(to^_, m_) — 5u(m_, m , )) = cos A5d; 
s = sin((5£)(m_|_, m_) — #D(m_, rn + )) — sin A8d- 



(5) 
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The phase difference ASu can be obtained from the 
sample of D mesons in a CP-eigcnstatc, decaying to 
Kgir + ir~. The Dalitz plot density of such decay is 

Pcp = |/d ± Id\ 2 = Pd +Pd ± 2 y / p D p D c (6) 

(the normalization is arbitrary). Decays of D mesons 
in CP eigenstate to Kgir + ir~ can be obtained in the 
process, e.g. e + e~ — > 0(3770) — ► DD, where the other 
(tag-side) D meson is reconstructed in the CP eigenstate, 
such as K + K~ or KgOj. 

Another possibility is to use a sample, where both D 
mesons (we denote them as D and D 1 ) from the ■0(3770) 
meson decay into the K a Tr + ir~ state [8j|. Since 0^(3770) 
is a vector, two D mesons are produced in a P-wave, 
and the wave function of the two mesons is antisymmet- 
ric. Then the four-dimensional density of two correlated 
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Dalit z plots is 

Pcorr(m%,m 2 _,m'*, m' 2 ) = \f D f' D ~ I'dId^ = 

I — : , , ( ? ) 

Pdp'd+PdPd ~ 2 Y PdPdp' d p' d {cc + ss), 

This decay is sensitive to both c and s for the price of 
having to deal with the four-dimensional phase space. 

In a real experiment, one measures scattered data 
rather than a probability density. Two options of deal- 
ing with real data are possible: a binned approach, or a 
scatter plot smoothing using nonparametric density esti- 
mation. The latter could be useful to reach the statistical 
sensitivity equivalent to the model-dependent case. How- 
ever, in this paper we show that using the appropriate 
binning this is also possible. 

III. BINNED ANALYSIS WITH Dcp DATA 

The binned approach was proposed by Giri et al. 
Assume that the Dalitz plot is divided into 2Af bins sym- 
metrically to the exchange m_ <-> m^_. The bins are de- 
noted by the index i ranging from —Af to Af (excluding 
0); the exchange m+ <-> m_ corresponds to the exchange 
i <-*• —i. Then the expected number of events in the bins 
of the Dalitz plot of D decay from B ± — ► DK ± is 

(Ni) = h B [K, + r%K-i + L\ 7vX~:.rr, + y Si )], (8) 

where Ki is the number of events in the bins in the Dalitz 
plot of the D° in a flavor eigenstate, h B is the normal- 
ization constant. Coefficients c, and Sj, which include 
the information about the cosine and sine of the phase 
difference, are given by 

/ \/pdPd cos( A6 D (m\,m 2 _))dV 

c = ; , (9) 

/ / PD dV j p D dV 

Si is defined similarly with cosine substituted by sine. 
Here T>i is the bin region, over which the integration is 
performed. Note that a = C-i, s, = — s_, and cf + sf < 1 
(the equality cf + sf = 1 being satisfied if the amplitude 
is constant across the bin). 

The coefficients Ki arc obtained precisely from a very 
large sample of D° decays in the flavor eigenstate, which 
is accessible at S-factories. The expected number of 
events in the Dalitz plot of Dcp decay equals to 

(Mi) = h C p[Ki + K_i + 2^/K~K~a\, (10) 

and thus can be used to obtain the coefficient q. As soon 
as the Ci and Si coefficients are known, one can obtain x 
and y values (hence, 03 and other related quantities) by 
a maximum likelihood fit using equation l[8]). 

Note that now the quantities of interest x and y (and 
consequently ^3) have two statistical errors: one due to 



a finite sample of B ± — > DK ± data, and due to Dcp — > 
Kgir + ir~ statistics. We will refer to these errors as B- 
statistical and -Dcf-statistical, respectively. 

Obtaining Si is a major problem in this analysis. If the 
binning is fine enough, so that both the phase difference 
and the amplitude remain constant across the area of 
each bin, expressions © reduce to Cj = cos(A#d) and 
Si = sin(A#D), so Si can be obtained as Sj = ±yl — cf. 
Using this equality if the amplitude varies will lead to the 
bias in the x, y fit result. Since Cj is obtained directly, and 
Si is overestimated by the absolute value, the bias will 
mainly affect y determination, resulting in lower absolute 
values of y. 

Our studies 0] show that the use of equality cf + s 2 = 1 
is satisfactory for the number of bins around 200 or more, 
which cannot be used with presently available Dqp data. 
It is therefore essential to find a relatively coarse binning 
(the number of bins being 10-20) which a) allows to ex- 
tract Si from Ci with low bias, and b) has the sensitiv- 
ity to the (f>3 phase comparable to the unbinned model- 
dependent case. 

Fortunately, both the a) and b) requirements appear to 
be equivalent. To determine the B-statistical sensitivity 
of a certain binning, let's define a quantity Q — a ratio 
of a statistical sensitivity to that in the unbinned case. 
Specifically, Q relates the number of standard deviations 
by which the number of events in bins is changed by 
varying parameters x and y, to the number of standard 
deviations if the Dalitz plot is divided into infinitely small 
regions (the unbinned case): 




(11) 

where f B = f D + {x + iy)f D , N t = \f B \ 2 dV. 

Since the precision of x and y weakly depends on the 
values of x and y 0], we can take for simplicity x = y = 0. 
In this case one can show that 

Q 2 u =y=0 = + sim y. n > ( 12 ) 

i I i 

Therefore, the binning which satisfies cf + sf = 1 (i.e. the 
absence of bias if Si is calculated as y/l — cf) also has the 
same sensitivity as the unbinned approach. The factor 
Q defined this way is not necessarily the best measure 
of the binning quality (the binning with higher Q can be 
insensitive to either x or y, which is impractical from the 
point of measuring ^3), but it allows an easy calculation 
and correctly reproduces the relative quality for a number 
of binnings we tried in our simulation. 

The choice of the optimal binning naturally depends 
on the D° model. In our studies we use the two-body 
amplitude obtained in the latest Belle 03 Dalitz analysis 

1- 
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FIG. 1: Divisions of the D° -> K° s n + ir- Dalitz plot. Uniform 
binning of ASd strong phase difference with M = 8 (left), and 
the binning obtained by variation of the latter to maximize 
the sensitivity factor Q (right). 

From the consideration above it is clear that a good 
approximation to the optimal binning is the one obtained 
from the uniform division of the strong phase difference 
ASd- In the half of the Dalitz plot rn\_ < m_ (i.e. the 
bin index i > 0) the bin T>i is defined by condition 

2n{i-l/2)/Af < A6 D (m 2 + ,m 2 _) < 2*(i + l/2)/AT, (13) 

in the remaining part (i < 0) the bins are defined sym- 
metrically. We will refer to this binning as A<5^-binning. 
As an example, such a binning with Af = 8 is shown 
in Fig. Q] (a). Although the phase difference variation 
across the bin is small by definition, the absolute value 
of the amplitude can vary significantly, so the condition 
cj + sf = 1 is not satisfied exactly. The values of a and 
St in this binning are shown in Fig. [2] with crosses. 

Figure Q] (b) shows the division with Af = 8 obtained 
by continuous variation of the Afe-binning to maximize 
the factor Q. The sensitivity factor Q increases to 0.89 
compared to 0.79 for A#£> -binning. 

We perform a toy MC simulation to study the statis- 
tical sensitivity of the different binning options. We use 
the amplitude from the Belle analysis [4j to generate de- 
cays of flavor D°, Dqp, and D from B ± -> DK ± decay 
to the K^tt + it~ final state according to the probability 
density given by ((T|), © and respectively. To ob- 
tain the ^-statistical error we use a large number of D° 
and Dcp decays, while the generated number of D de- 
cays from the B ± — > DK^ process ranges from 100 to 
100000. For each number of B decay events, 100 samples 
are generated, and the statistical errors of x and y are 
obtained from the spread of the fitted values. A study 
of the error due to Dcp statistics is performed similarly, 
with a large number of B decays, and the statistics of 
Dcp decays varied. Both errors are checked to satisfy 
the square root scaling. 

The binning options used are A^-binning with Af = 8 
and Af = 20, as well as "optimal" binnings with maxi- 
mized Q obtained from these two with a smooth variation 
of the bin shape. Note that the "optimal" binning with 
Af = 20 offers the ^-statistical sensitivity only 4% worse 
than an unbinned technique. For comparison, we use the 
binnings with the uniform division into rectangular bins 



(with Af = 8 and Af = 19 in the allowed phase space, the 
ones which are denoted as 3x3 and 5x5 in ||). 

The B- and Ucp-statistical precision of different bin- 
ning options, recalculated to 1000 events of both B and 
Dcp samples, as well as their calculated values of the 
factor Q, are shown in Table [U In the present study we 
use the errors of parameters x and y rather than <f>^ as 
a measure of the statistical power since they are nearly 
independent of the actual values of 03, strong phase S 
and amplitude ratio tq. The error of 4>3 can be obtained 
from these numbers given the value of tb- The factor 

Q reproduces the ratio of the values ^j\ja\ + 1/c^ for 

the binned and unbinned approaches with the precision 
of 1-2%. While the binning with maximized Q offers 
better £?-statistical sensitivity, the best Dcp-statistical 
precision of the options we have studied is reached for the 
A(^-binning. However, for the expected amount of ex- 
perimental data of B and Dcp decays the ^-statistical 
error dominates, therefore, slightly worse precision due 
to Dcp statistics does not affect significantly the total 
precision. 

Using AfJo-binning, the following combination of the 
binned and unbinned approaches is possible, which al- 
lows to reach B-statistical precision equivalent to the 
unbinned case. Assume the number of A<5o-bins is large 
enough, so cos A<5d and sin A5d remain almost constant 
across the bin area. At the level of current precision this 
is reached already for a number of bins as small as 10-20. 
Then 

Cj = cosA^d — , (14) 

/ / PD dV J p D dV 

where the integrals can be calculated from the flavor D 
data sample. Therefore, it is possible to obtain cos A^d 
from a (and consequently sin ASd), and use them in the 
expression for the probability density ^ to perform the 
unbinned £>-data fit thus obtaining the best possible B- 
statistical precision. 

We have considered the choice of the optimal binning 
only from the point of statistical power. However, the 
conditions to satisfy low model dependence are quite dif- 
ferent. Since the bins in the binning options we have 
considered arc sufficiently large, the requirement that the 
phase does not change over the bin area is a strong model 
assumption. We have performed toy MC simulation to 
study the model dependence. While the binning was kept 
the same as in the statistical power study (based on the 
phase difference from the default D° amplitude), the am- 
plitude used to generate D°, D C p and B ± -> DA" ± de- 
cays was altered in the same way as in the Belle study 
of the model-dependence in the unbinned analysis [J]. 
As a result, the same bias of A03 ~ 10° is observed as 
in unbinned analysis. We remind that the cause of this 
bias is a fixed relation between the Cj and Sj. Therefore, 
proposed binning options, although providing good sta- 
tistical precision, are not flexible enough to provide also 
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TABLE I: Statistical precision of (x, y) determination using different binnings and with an unbinned approach. The errors 
correspond to 1000 events in both the B and Dcp ((K°irir) 2 ) samples. 
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a low model dependence. To minimizie the model depen- 
dence, the bin size should be kept as small as possible, 
therefore, uniform binning is more preferred. 

In a real analysis, one can control the model error by 
testing if the amplitude used to define binning is com- 
patible with the observed Dcp data. This can be done, 
e.g., by dividing each bin and comparing calculated val- 
ues of Ci in its parts, or by comparing the expected and 
observed numbers of events in each bin. 

We conclude that the method of 3 determination 
using only Dqp data is only asymptotically model- 
independent, since for any finite bin size the calculation 
of Si is done using model assumptions of the ASd varia- 
tions across the bin. Increasing the Dcp data set, how- 
ever, allows to apply a finer binning and therefore reduce 
the model error due to the variation of the phase differ- 
ence. 



IV. BINNED ANALYSIS WITH CORRELATED 
D° -> K°irn DATA 

The use of the "0(3770) decays where both neutral D 
mesons decay to the ATg7r + 7r _ state allows to significantly 
increase the amount of data useful to extract phase infor- 
mation in D° decay. It is also possible to detect events 
of 0(3770) (K s n + n-) D (K L 7r+7r-) D , where K° L is not 
reconstructed, and its momentum is obtained from kine- 
matic constraints. The number of these events is approx- 
imately twice that of (Kgirir) 2 due to combinatorics. We 
will refer to both of these processes as (K°irir) 2 . 

In the case of a binned analysis, the number of events 
in the region of the (Kgirir) 2 phase space is 

{M)ij = HcorrlKiK-j + K^Kj- 

' _ (15) 

1\ ,[(■,(■, + SiSj)]. 

In the case of the (Kgir + Tr~)(Klir + Tr~~) final state, the 
interference term changes its sign. Here two indices corre- 
spond to two D mesons from ^(3770) decay. It is logical 
to use the same binning as in the case of Dcp statis- 
tics to improve the precision of the determination of 
coefficients, and to obtain Sj from data without model 
assumptions, contrary to Dcp case. The obvious advan- 
tage of this approach is its being unbiased for any finite 



(K°irir) 2 statistics (not asymptotically as in the case of 
Dcp data). 

Note that in contrast to Dcp analysis, where the sign 
of Si in each bin was undetermined and has to be fixed 
using model assumptions, (Kgir + ir~) analysis has only a 
four- fold ambiguity: change of the sign of all Ci or all Si. 
In combination with Dcp analysis, where the sign of c,j 
is fixed, this ambiguity reduces to only two-fold. One of 
the two solutions can be chosen based on a weak model 
assumption (incorrect Si sign corresponds to complex- 
conjugated D decay amplitude, which violates causality 
requirement when parameterized with the Breit-Wigncr 
amplitudes). 

Coefficients Ci, Sj can be obtained by minimizing the 
negative logarithmic likelihood function 

- 2 log£ = log P(M i:i , (M)ij), (16) 

id 

where P(M, (M)) is the Poisson probability to get M 
events with the expected number of (M) events. 

The number of bins in the 4-dimcnsional phase space 
is 4Af 2 rather than 2Af in the Dqp case. Since the ex- 
pected number of events in correlated KgTrir data is of 
the same order as for Dcp, the bins will be much less 
populated. This, however, does not affect the precision 
of Ci, Si determination since each of the free parameters 
is constrained by many bins. 

The toy MC simulation was performed to study the 
procedure described above. Using the amplitude from 
the Belle analysis, we generate a large number of D° — > 
KgTr + TT~ decays and several sets of (K^iiit) 2 decays (ac- 
cording to the probability density given by ([7])). We use 
the same binning options as in Dcp study with Af = 8. 
The negative logarithmic likelihood (fT6|) is then mini- 
mized with Ci, Si and h cmr as free parameters. We con- 
strain |ci,Si| < 1 in the fit to improve the convergence. 
The coefficients Ci, Si are then used in the fit to B de- 
cay data to obtain the x,y error due to (K°irir) 2 decay 
statistics. The number of (KgTnr) 2 decays ranges from 
1000 to 10000. The obtained error shows a square root 
scaling. The best (if°7T7r)-statistical error is obtained 
for A<5D-binning and recalculated to 1000 events yields 
a x = 0.0050, (J y = 0.0095, which is only slightly worse 
than the error obtained with the same amount of Dqp 
data (see Table Q] for comparison) . We also check that 
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FIG. 2: Coefficients a, Si calculated from 2000 generated 
(K°mv) 2 events (points). A&> binning with Af = 8, 100 toy 
MC samples are used. Different colors correspond to different 
bins. Crosses denote the true Ci, s* values in each bin. 

changing the model used to define the binning does not 
lead to the systematic bias (although it does decrease 
the statistical precision). Figure [5] demonstrates the pre- 
cision of the determination of c,, Sj coefficients in our 
toy MC study and the absence of the systematic bias for 
both Ci and s,. 

Since the number of (K°inr) 2 decays in 0(3770) data 
is approximately twice larger than the number of Dcp 
decays, the statistical errors due to ^(3770) data for 
the two approaches are nearly equal. The same binning 
can be used in both approaches, therefore improving the 
accuracy of c% determination. The approach based on 
(K°tttt) 2 data allows to extract both a and s, without 
additional model uncertainties, so it can be used to check 
the validity of the constraint cf + sf = 1 and therefore to 
test the sensitivity of the particular binning. 

V. CONCLUSION 

We have studied the model-independent approach to 
03 measurement using B ± — ► DK ± decays with neutral 



D decaying to K^ti + -k" . The analysis of ^(3770) — > 
DD data allows to extract the information about the 
strong phase in D — > Kgir + ir~ decay that is fixed by 
model assumptions in a model-dependent technique. We 
specially consider the case with a limited ?/>(3770) — > DD 
data sample which will be available from CLEO-c in the 
nearest future. 

In the binned analysis, we propose a way to obtain 
the binning that offers an optimal statistical precision 
(close to the precision of an unbinned approach). Two 
different strategies of the binned analysis are considered: 
using Dcp — ► A'^tt+tt - data sample, and using decays 
of ^(3770) to (K°Tr + TT-) D (K°Tr + n-) D . The strategy us- 
ing Dcp decays alone cannot offer a completely model- 
independent measurement: it provides only the informa- 
tion about Ci coefficients, while Si for low Dcp statistics 
has to be fixed using model assumptions. However, as 
the Dcp data sample increases, model- independence can 
be reached by reducing the bin size. The strategy us- 
ing the V(3770) — > (K Tr + TT~)ri(K ir + TT~)D sample, in 
contrast, allows to obtain both Ci and with an accu- 
racy comparable to Dcp approach. Both strategies can 
use the same binning of the D° — > Kgir + TT~ Dalitz plot 
and therefore can be used in combination to improve the 
accuracy due to ^(3770) statistics. 

The expected sensitivity is obtained based on the 
D° decay model from Belle analysis. For the CLEO- 
c statistics of 750 pb" 1 (1000 D C p events and 2000 
(K°irir) 2 events) the expected errors of parameters x 
and y due to '0(3770) statistics are a x = 0.003 and 
(jy = 0.007. For tb = 0.1 it gives the 03 precision 
(703 = max((T I , <jy) / (\J~2rB) — 3°, which is far below the 
expected error due to present-day B data sample. 

In our study, we did not consider the experimental sys- 
tematic uncertainties e.g. due to imperfect knowledge of 
the detection efficiency or background composition. We 
believe these issues can be addressed in a similar manner 
as in already completed model-dependent analyses. 
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